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Abstract 

We show that the hnearized equations of the incompressible elastic medium admit a 'Maxwell form' in which the shear 
component of the stress vector plays the role of the electric field, and the vorticity plays the role of the magnetic field. 
Conversely, the set of dynamic Maxwell equations are strict mathematical corollaries from the governing equations of 
the incompressible elastic medium. This suggests that the nature of 'electromagnetic field' may actually be related to 
an elastic continuous medium. The analogy is complete if the medium is assumed to behave as fluid in shear motions, 
while it may still behave as elastic solid under compressional motions. Then the governing equations of the elastic 
fluid are re-derived in the Eulerian frame by replacing the partial time derivatives by the properly invariant (frame 
indifferent) time rates. The 'Maxwell from' of the frame indifferent formulation gives the frame indifferent system 
tiiat is to replace the Maxwell system. This new system includes terms already present in the classical Maxwell 
equations, alongside terms that are the progenitors of the Biot-Savart, Oersted-Ampere's, and Lorentz-force laws. 
Thus, a frame indifferent (truly covariant) formulation of electromagnetism is achieved from a single postulate: the 
electromagnetic field is a kind of elastic (partly liquid, partly solid) continuum. 

Keywords: Frame Indifference, Maxwell's Elastic Fluid, Maxwell's Equations of Electrodynamics, Lorentz Force, 
Biot-Savart Law 



1. Introduction 

Classically, wave motion in elastic media is a phe- 
nomenon associated with either/both the transverse or 
longitudinal vibrations; if there is a wave, then some- 
thing material should be waving. This notion led 19th 
century scientists to introduce the concept of the lu- 
miniferous medium (field, aether, etc.). The first at- 
tempt to explain the propagation of light as a field phe- 
no menon was by C auchy circa 1827 (see the account 
in (IWhittakerl ll989h X who postulated the existence of 
an elastic continuum through which light propagates as 
a shear wave. Subsequently came the contributions of 
Faraday and Ampere, which eventually led to the for- 
mulation of the electromagnetic model as known today. 
The crucial advance in understanding the phenomena of 
electr omagneti s m ph enomena was achieved, however, 
when lMaxwelll (1 1 865h added the 'displacement current' , 



to Ampere's law. This term was similar to the time- 
derivative of the stress in his constitutive relation for 
elastic gases (' Maxwell . 1867 ). Since the electric field 
vector is a clear analogue of the stress vector in con- 
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tinuum mechanics, one can say that Maxwell postu- 
lated an elastic constitutive relation by adding the dis- 
placement current to Ampere's law. Indeed, the new 
term transformed the system of equations (akeady es- 
tablished at that time in electrostatics) into a hyperbolic 
system with a characteristic speed of shear wave propa- 
gation. Maxwell and Hertz identified the characteristic 
speed with the speed of light. 

However, soon after Maxwell formulated his equa- 
tions, it was discovered that his model was not invari- 
ant with re spect to trans lational motion of the coordi- 
nate frame. iHertzld 19001) realized that the cause of non- 
invariance was the use of partial time derivatives. He 
proposed to use the convective time derivative in lieu 
of the former The Maxwell-Hertz equations (called 
also 'progressive wave equations') are clearly the cor- 
rect model for electromagnetic phenomena in moving 
bodies. But the primordial question is whether the 
progressive-wave equations can be construed to hold 
also in vacuo. The answer is obviously in the affirma- 
tive if one accepts that what is currently called 'physical 
vacuum' must be regarded as a material continuum. 

Voigt (see lErnst and Hsul ( l200ll) ). lLarmod(ll897l ), and 
independently Lorentz, spotted the fact that the wave 
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equation can be made invariant if, in the moving frame, 
the time variable is changed together with the spatial 
variables. Nowadays, this is known as the Lorentz 
transformation. The success of the latter stems from 
the fact that it tacitly restores some parts of the con- 
vective derivative, i.e. it emulates the material invari- 
ance for non-deformab le frames in rectilinear motion 
dChristov , I 2006ll2009bl) . Although researchers usually 
speak about Lorentz covariance as a general covariance, 
it has to be pointed out that the Lorentz transformation 
has no meaning for accelerating frames (nor for gener- 
ally deforming frames). Hence the search for the truly 
covariant formulation should continue. 

The present paper summarizes the efforts from the 
last decade and a half towards identifying the me- 
chanical construct behind the phenomenon of what is 
termed 'electromagnetic field'. We show here that 
the Maxwell equations and the laws of electromagnet- 
ism (Biot-Savart, Oersted-Ampere, and Lorentz-force 
laws), all can be derived from the governing equations 
of Maxwell's elastic fluid, that includes relaxation of 
the stress. In doing so, we also present a concise 
frame indifferent formulation of the Maxwell elastic 
fluid model. 



2. Maxwell Form of the Equations of Incompress- 
ible Elastic Medium 

The linearized governing equations of elastic medium 
are valid only for infinitesimal deformations, when the 
referential description X and spatial configuration x co- 
incide. The linear constitutive relationship for an elastic 
body relates the stress tensor & (we reserve the notation 
cr for the deviatoric stress tensor) to the deformation ten- 
sor e via the generalized Hooke's law 

a-^ Atr(e) + 27]e, e ^[Vh -H (Vh)T], (1) 

where tr(e) is the trace of the deformation tensor, and 
the superscript j denotes the transpose. The above con- 
stitut ive law yield s the so-called Navier equations (see, 
e.g., dSegel , 1987 , pg. 1 17) for the displacement vector 
uix; t): 

tArr = + '7)V(V ■ h) + T]V\ 



dfi 



^{A + 277)V(V ■u)-j]VxVxu, (2) 



where A and // are the Lame coefficients. Here one can 
use the 'nabla' operator V because the equations are 
written in the current configuration. 

The speeds of propagation of shear and compres- 
sional disturbances, are given respectively by 



2rj + A 



lrj + A 



(3) 

ci 



where the ratio 6 is introduced for convenience. 

In a compressible elastic medium, both the shear 
and the dilational/compressional waves should be ob- 
servable. Since the groundbreaking works of Young 
and Fresnel, it is well established that electromagnetic 
waves (e.g., light) are a purely transverse (shear) phe- 
nomenon. This observation requires us to reduce the 
complexity of the model and to find a way to elimi- 
nate the dilational modulus A. Cauchy assumed that 
A - and e nded up wit h the th eory of so-called 'volatile 
aether' (see IWhittaken (Il989h ). Upon a closer exami- 
nation, we found that such an approach cannot explain 
Maxwell's equations. 

Let us now assume that the dilational waves are not 
observable because the other extreme situation is at 
hand: A^ rf which is equivalent to 5 <s 1 or (5 ' » 1. 
It is convenient to rewrite Eq. ^ in terms of the speed 
of light, c, and parameter 6, namely 



,d u 



5(c"^— -H V X V X m) = V(V ■ m), 

^ or ' 



(4) 



and to expand the speed of light c, displacement u, and 
velocity v :- du/dt into asymptotic power series with 
respect to 6, namely 

c - Co + Sci + ■ ■ ■ , 

u = uq + 6ui + ■ ■ ■ , (5) 

V = vo + (5vi + ■ ■ ■ . 

Introducing Q into dUi and combining the terms with 
like powers we obtain for the first two terms 



V(V.Ho) = 0. 



(6a) 



_2 d^uo 



def 



1 



-H V X V X Mo = V(V • Ml) = -Vp , (6b) 



where p is introduced for convenience and is the coef- 
ficient of the spherical part of the internal stresses. The 
variable p has dimension of jic^ and plays the same role 
as the pressure in an incompressible medium: p is an 
implicit function in Eq. ( l6b] i that provides the necessary 
degree of freedom to enforce the satisfaction of the 'in- 
compressibility' condition, Eq.(l6a]l. The latter can also 
be rewritten as 



V • M() = const. 



V ■ Vo = 0, 



(7) 



which requires that the velocity field be solenoidal 
within the zeroth-order of approximation of the small 
parameter 6. From now on, the subscript '0' will be 
omitted form the variables without fear of confusion. 
Now, Eq. (l6bl l can be rewritten as 



def 



S = V ■ cr = -77V X V X H, (8) 
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where s is the tangential part of the stress vector in the 
elastic continuum. The normal part of the stress vector 
is given by the pressure gradient Vp. 

Let us now introduce the following notations 



def 



-s = 77 V X (V X m), 

def _ 7_ 

B - i^V XV - T]c V X V, 



(9a) 
(9b) 



and call them the 'electric field' and 'magnetic induc- 
tion' vectors. Thus the electric field is the negative tan- 
gential stress vector, while magnetic field is the 'dy- 
namic' vorticity. Naturally, //:= Vxv is the 'kinematic' 
vorticity. In the virtue of the above definition, one has 
the following equation for the magnetic field 



V ■ B = 0. 



(10) 



Now taking the curl of Eq. dHJ and acknowledging the 
definitions given in Eqs. (|9]l, we arrive at Faraday's law 



(11) 



„ ^ dB 
VxE = . 

dt 



On the other hand, taking the time derivative of Eq. (I9al i. 
we obtain the second of the dynamic Maxwell's equa- 
tions 



dE , 
— = 77 V X (V x v) = c'^V X B. 
of 



(12) 



The fact that the governing equations of any elastic 
medium in the linear limit admit a Maxwell form also 
can be considered as an indication that the electromag- 
netic field is in itself an elastic medium. In what follows, 
we shall call the mechanical object equivalent mathe- 
matically to the electromagnetic field the metacontin- 
uum, to distinguish it from continuous media in techni- 
cal applications, such as fluids and elastic solids. Note 
that the inverse of the elastic shear coefficient plays the 
role of the electric permittivity in vacuo, while the den- 
sity of the metacontinuum acts as the 'magnetic perme- 
ability.' 

The results of this section unequivocally show that 
the 'field' described by Maxwell's equations is equiv- 
alent to an elastic material. To best of author's knowl- 
edge, the connection of Maxwell's equation to the equa- 
tions governing elastic media was first established in 
(Christov, 1995, 1996, 2001)Q The common theme of 
these earlier papers is that the metacontinuum is consid- 
ered as an elastic solid. In such a model, no stationary 
magnetic fields can exist, since no steady velocities are 
possible for a solid continuum without discontinuities. 



' Similar deri vations were proposed, presumably independently, in 
tomitrievll200l . 



3. The Electromagnetic Field: Liquid or Solid? 

In this section we outline the next decisive step in de- 
veloping the model: we consider an elastic liquid in lieu 
of an elastic solid. This means that for shear deforma- 
tions, the metacontinuum must be an elastic fluid for 
which the time derivative of the deviatoric stress tensor 
is related to the deviatoric rate of deformation tensor x 
via the relation 



da- 



X = i(Vv + Vv^) - V ■ 



(13) 



This constituitve relationship (rheology) can be rewrit- 
ten for the negative stress vector E - -V ■ cr and defor- 
mation vector, d -V ■ X {- -V x V x v), since both of 
these vectors are the divergences of the respective ten- 
sors involved in the elastic rheology, Eq. (fTjt . Then 



dE dE 
= ^ V X V X V, or 



dt 



dt 



T]VxVxv, (14) 



where ^ is called 'elastic viscosity' by iJoseph et al 
(1986), T is the relaxation time of the stresses, and the 
apparent elastic shear modulus is given by 77 = ^t"'. 
Note that the above elastic-liquid rheology concerns just 
the shear deformations. For compressional/dilational 
motions, the metacontinuum can still behave as a virtu- 
ally incompressible solid. A more general formulation 
of the shear part of the constitutive relation would be as 
in the viscoelastic liquid 

dE dE 
T — +ktE - (VxVxv, or — +kE - rjVxVxv, (15) 
dt dt 

where k can be called the 'cond uctivity' of the vis- 
coelastic liquid. Note that in (Jose ph et al.L 1986), the 
conductivity is set to unity, which precludes treating 
purely elastic (non-viscous) liquids. Here, we prefer to 
keep the flexibility offered by the presence of the coef- 
ficient K. Setting the appropriate terminology is an un- 
easy task because, when k 0, then ^ does indeed have 
a meaning of a viscosity coefficient, while for k - 0, 
it loses its independent meaning and enters the picture 
through the coefficient of apparent shear elasticity 77. 

The case k Q gives a viscoelastic rheology, but 
it does not lead to a model governed by the Navier- 
Stokes equation^ with additional elasticity, because no 
retardation term (time derivative of the deformation ten- 
sor/vector) is present. In this sense, adding the con- 
ductivity does not introduce dispersive dissipation, but 
rather a linear attenuation parameterized by the conduc- 
tivity coefficient. For the effects connected w ith the at- 
tenuation/conductivity we refer the reader to ( Christovl 
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2007h and jHarmuthl Il986h . The constitutive relation- 
ship given in Eq. ( fT5b can be interpreted as Ohm's law 
for vacuo. Althoug h, thi s stipu lation is made in main- 
stream texts (e.g., jjoosl 1986 )). the cause of Ohm's 
law in matter is not necessarily the intrinsic resistance 
of the metacontinuum: it is the result of the thermal 
fluctuations of the atoms that obstruct the free passage 
of charges through a conductor Clearly, a more in- 
depth argument is needed to justify having Ohm's law 
in vacuo, which goes beyond the scope of the present 
paper 

It should be pointed out here that Eq. (fTST i is con- 
cerned merely with the rheology for the shear motions. 
At this point it is not of importance if for dilational / 
compressional motions the metacontinuum is solid or 
liquid, provided that the dilational coefficient is much 
larger than the shear coefficient. Then, 



tr(,T) 



|(y + 20V-v 
(/1-h 277)V-h 



liquid, 
solid. 



(16) 



The intuitive argument is that if the metacontinuum is a 
liquid with respect to the dilational motions, it may lose 
its integrity during the motion. Since there is no infor- 
mation on the electromagnetic field being 'ruptured,' we 
are guided by the above argument and assume that the 
solid rheology applies to the dilational/compressional 
motions. Note that in the above notations, v is the dila- 
tional viscosity coefficient, and ^ is the shear viscosity 
coefficient defined in Eq. (fTSlli . 

The system (fTSl ). (fTSI l? specify the combined consti- 
tutive relation for the metacontinuum (electromagnetic 
field). The closest analogy to an ubiquitous contin- 
uous medium is that of jelly or pine pitch. If com- 
pression/dilation waves are sent through the metacon- 
tinuum, it behaves as a elastic solid with very large dila- 
tional modulus, while if a shear deformation is applied, 
it flows as a incompressible liquid. 

4. Euler Variables: the Frame Indifferent Model of 
Incompressible Maxwell Fluid 

In terms of the velocity vector, the Cauchy balance 
can be rewritten as follows: 

dv 

■ Vv = divo- = -V/? (17) 

at 

where the left-hand side is the material (convective) 
derivative of the velocity vector v in the current con- 
figuration (called 'convective' or 'total' derivative). Re- 
member that in the referential description, it is just the 
partial time derivative. Note also that for an incompress- 
ible metacontinuum, the density is the same constant in 



both the referential and spatial descriptions and we de- 
note it by yU. The concept of frame indifference (gen- 
eral covariance of the system) requires that the partial 
derivative of the stress variable (in our case the stress 
vector -E) in Eq. ( fTSl ) is replaced by the appropriate 
invariant rate. Since i t is a v ector density, see the argu- 
ment by Schrodinger ( 1950l) . the rate has to leave the 
integral of the s tress vector invariant. It is argued in 
dChristov , l 2006h that the pertinent invariant rate is the 
so-called Oldroyd upper-convected derivative, namely 



m def dE 



^1|f£^+v.V£-£.Vv + (V.v)£. 
at at 



(18) 



Here we come to one of the most crucial assumptions 
of the present work, namely, the way the constitutive 
relation has to be written when a relaxation of the stress 
is present. It is usually assumed that the invariant rate 
to be added to the constitutive law of viscous liquids 
should be of the stress tensor, i.e. the upper-convected 
Oldroyd derivative (.Oldrovd. (.1949) : . Bird et al..a987,) ): 

9cr dcr 

-— :— hv- V(T-cr grad v-(gradv)' cr -Her divv. (19) 

of ot 

The problem with this constitutive conjecture is that 
the tensor cr does not play an independent role in the 
Cauchy balance equation. Rather, the deviatoric stress 
vector -E := diver appears there. Then, does Eq. (fT9] l 
ensure that the time rate of div cr is invariant? To find 
out, we take the operation div of Eq. (fT9] l, arriving at 

,. do- 5div cr 
div(-5-) = — + V • V(div cr) 
at of 

- (grad v)^(div cr) + (div cr) div v - ( VVv)cr 
3 div cr 



■ (VVv)cr, (20) 

Of 

which differs from the invariant time rate of div cr by 
the term (- VVv)cr, the latter being the contraction of the 
third rank tensor of the repeated gradient (the Hessian) 
of V and the second rank stress tensor cr. In order to 
establish which constitutive relation is correct (the one 
involving the stress tensor or stress vector), one has to 
have data for flows in which the Hessian of the velocity 
vector field is measured. While for the case of elastic 
liquids it is still possible to devise such an experiment, 
electromagnetism clearly indicates that the constitutive 
relation at play is the one involving the stress vector 

Here we propose an alternative formulation for mod- 
els involving stress relaxation by replacing the partial 
time derivative of the stress vector with the invariant rate 

dE 

— +vVE-E-'Vv + E(V-v) + KE^nVx(Vxv), (21) 
ot 
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which was proposed in dChristov , 2006 ), and then suc- 
cessfully applied to the generalizati on of the Maxwe ll- 
Cattaneo model of second sound in dChristovi l2009al) . 

Apart from providing insight into the possible consti- 
tution of the electromagnetic field, the above model has 
practical significance for the theory of Maxwell elastic 
fluids. Note that we add the condition of incompress- 
ibility (recall Eq. Q with the subscript "0" omitted): 



V • V = 0. 



(22) 



Collectively, we can term the system Eqs. ( fTTl i. 
(ISTT l and ( |22] | as the 'Frame Indifferent Incompressible 
Maxwell Fluid Model' . 

5. Eliminating the Stress Vector 

The alternative formulation of the constitutive rela- 
tion based on the stress vector concept proposed here 
has a very important consequence: it allows the stress 
vector -E to be eliminated between Eqs. (l2Tl i and (fTTI l 
(when V • V = 0), to obtain a system that does not contain 
the stress variable : 

d^v dv (3V» 
^ — _ + 2^'i' ■ ^ 1- (wwxiv — - 



dt 



5f 



+ T/V^v, (23) 



where the following notations are used 



gV/7 



(vv)VVv:=^.^^2.^^y,y,V,V,v. 



V ■ V(V/?) - dp) ■ Vv, (24) 
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(25) 



Note that 3/9? is the upper convected derivative of a 
vector defined in Eq. ( fTsT l. 

Eq. ( |23] ) contains the implicit function p that has to 
ensure the satisfaction of the incompressibility condi- 
tion Eq. ( l22l l. The advantage here is that, the only un- 
known functions are v and p, and the hyperbolicity of 
the model is now easily seen. 

6. Frame Indifferent Electromagnetodynamics 

According to the principal of frame indifference 
dTruesdelli [19651), the laws of physics (including the 
laws of continuum physics) must have the same form in 
any reference frame (coordinate system). Unlike what 
is called 'Lorentz Covariance', the laws in the refer- 
ential description are frame indifferent, i.e. they are 
truly covariant. However, experimental measurements 
are always connected with a current frame in the geo- 
metric space. This means that an observational frame 
cannot detect the material variables (the referential de- 
scription), but rather can merely measure their counter- 
parts in the current (geometric) frame. This is a typi- 
cal situation in mechanics of continuous media where 



the reference configuration is often not related to any 
measurable frame. For this reason we need to reformu- 
late the model from Section|2]in the current description 
making use of Euler variables. This is the objective of 
the present section. 

The Cauchy balance, Eq. ( [TtI i. can be rewritten in the 
so-called 'Gromeka-Lamb form' (Sedov, 1981): 



fi—-l^vx(Wx v) 
of 



-V(p - E, 



p + y, (26) 



Now, taking the curl of Eq. 
nitions Eq. (|9]l, we get: 



, and using our defi- 



Vx[£-vxB] = 

of 



(27) 



which is Faraday's law with an additional term repre- 
senting the force exerted by a moving magnetic field on 
each point of the medium. It is induced by the con- 
vective part of the acceleration at that point. The reac- 
tion to this force is the force acting on a moving point 
(charge), v x B, known as the 'Lorentz force'. In other 
words, the material invariant version of Faraday's law 
presented here automatically accounts for the physical 
mechanism that causes the Lorentz force. The latter is 
nothing more than the inertial force given by the con- 
vective part of the total derivative. This is a very impor- 
tant result because it tells us that the Lorentz force is not 
an additional, empirically observed force that has to be 
grafted onto Maxwell's equations, but rather, it is con- 
nected to the material time derivative, specifical, to its 
convective part. Under the incompressibility condition, 
Eq. d27b can be recast as 

dB 

— + V ■ VB - B ■ Vv = -V X E, (28) 
dt 

which we can call the 'Hertz form' of the Faraday- 
Lorentz law. Note the presence of the third term on the 
left-hand-si de. It is not in Hertz's f ormulation, nor does 
it appeal- in dChristovi I2OO6I l2007b . Evidendy, Eq. ^ 
does not give any special advantage over Eq. d27] i. but 
the form of Eq. d28] l shows that one cannot just add 
the convective part of the derivative to Faraday's law 
to make it frame indifferent. This follows from the fact 
that the magnetic field is not a primary variable (primary 
variables in fluids are the velocity components and the 
pressure). Rather, it is proportional to the curl of the 
velocity vector, v. 

By using the vector identity (see, e.g., 
dBorisenko and Tarapov , Il979[ pg. 180) 



Vx(vx£) = £-Vv-v-V£-Hv(V-£)-£(V-v), (29) 

Eq. d2TI) yields the following generalization of the sec- 
ond of the dynamical Maxwell's equation: 

^ - V X (v X £) +kE = -v(V^) -I- c-WxB , (30) 
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where the definition given in Eq. ( l9b] l is already ac- 
knowledged. 

Now, Eqs. (|27ll, (HE]), ([TOll form a system which 
can be termed the Frame Indifferent Electromagneto- 
dynamics (FIEM). This system generalizes Hertz's pro- 
gram from 1890 and rigorously fulfills the requirements 
for 'General Covariance,' because FIEM is frame- 
indifferent; it is invariant when changing to another co- 
ordinate frame that can accelerate and can even deform. 
A very simple, limiting case of frame indifference is that 
of Galilean invariance. The vector of absolute velocity 
V is the primary variable, but the absolute velocity in the 
referential description cannot be measured (only rela- 
tive velocities can be observed). In principle it can be 
restored from the magnetic field Eq. (|9bt . provided a 
boundary condition is known. 

Clearly, in the limit of small velocities the convective 
and convected terms can be neglected and in the limit 
one obtains Maxwell's system. 

7. Interpretation of the Different Terms in the Con- 
vected Time Rate 

A remarkable feature of Eq. (l30t is that it incorpo- 
rates terms (with under-braces) that, collectively, form 
the progenitor of the Biot-Savart law. Indeed, set- 
ting aside the possible singularities connected to point 
charges, then we can neglect the term with V ■ E. Let us 
also consider the case of stationary electric field, E, - 0, 
and no resistance k - Q. Then, considering a surface D 
inside a closed contour C, we can integrate Eq. (|30] | over 
C and use the Stokes theorem to obtain 



/X 



[vxE + B]dis^O, 



(31) 



where ds is an areal element on the surface D. Since the 
surface D is arbitrary, then the only possibility is that 
the integrand must be equal to zero. Thus one arrives 
at the form of the Biot-Sava rt law as it is s tipulated in 
relativistic electrodynamics (iGriffiths[|l98U §12.3), 



B 



1 



-vxE 



(32) 



where v is the velocity of the point of the field at which 
B is measured. 

The underlined terms in Eq. (|30] | give that in vacuo 



VxB^^E, (33) 



which is discussed in (iJoos . Ch.l5) in dimen- 

sional form. This can be interpreted as a combination 
of Ohm's and Ampere's laws of electromagnetism in 
vacuo. 



The terms with hats over them in Eq. (|30) give as a 
corollary the Ampere law if the following definition of 
a charge in vacuo is introduced: 



def 

p = V 



(34) 



In particular, the 'chargedness', p, of the displace- 
ment/velocity field is defined as the divergence of the 
electric field at the point. In order not to confuse this 
property of the field in vacuo with the localized pattern, 
called electron or a proton, we call the above defined 
function p the metacharge. For the latter, a continuity 
equation is readily derived upon applying the operation 
div to Eq. ( l30l ). namely 



dp 
dt 



+ V ■ (pv) + Kp -Q, 



(35) 



which, for the case k - 0, \s the standard continu- 
ity equation for charge. Here, is to be mentioned that 
Eq. ( |35] | was derived by Chri stov (2006) directly from 
the Oldroyd form, Eq. (fTsT l. but the derivation here is 
much more straightforward because of the application 
of the identity Eq. (|29] |. We shall refer to Eq. (l35t as the 
continuity equation for the metacharge. 

In terms of the above introduced metachargJ^ we ob- 
tain 



c^V X B = pv = /, 



(36) 



which can be called Ampere's law in vacuo. 

The important conclusion form the frame-indifferent 
formulation of the displacement current is that simi- 
larly to the Lorentz-force law, the convective/convected 
terms are related to phenomena that are embod- 
ied in Ampere's and Biot-Savart's laws, thus uni- 
fying them with Maxwell's electrodynamics. All 
three electromagnetic-force laws (called alternatively 
the Taws of motional electromagnetism') are manifes- 
tations of the inertial forces in the metacontinuum. 



Whittaked(ll989l pp. 85-88) put forward the idea that 



both of these laws may actually follow from a single 
law, similar to what is presented in Eq. (l30t . A de- 
bate is still ongoing in the htera ture about wether these 
two laws are identical (see, e.g., ("ChristodouUdes',' 19871 
Jollv, 1985)) or independent (see, e.g., Graneau ( 1985)). 
Our results seem to favor Whittaker's original idea that 
both laws have to be interwowen in the correct formu- 
lation. In our work they are merely the corollaries form 
the inertial terms embodied in the convected derivative. 



^^Note that the ubiquitous notion of charge can be explained away 
by a localized shear pattern in the metacontinuum, and then the usual 
charge is the integral of the metacharge over the span of the localized 
pattern fsee iChristovi t2007i) ). 
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The frame indifferent model of the electromagnetic 
field (called here the metacontinuum), succeeds in uni- 
fying, in a single nexus, all known phenomena of 
electromagnetism: Faraday's law, displacement-current 
law, Lorentz-force law, Oersted-Ampere's law, and 
Biot-Savart law. It is a significant step forward from 
Maxwell's model, in which only the first two were ex- 
plained by the equations themselves, and the latter three 
appeared as additional empirically observed relations 
between the main characteristics of the field. Since 
in our model these new terms are valid in vacuo, they 
are clearly the progenitors of the related phenomena in 
moving media. 

8. Effect of Compressibility (Dark Energy?) 

As shown above, the model that yields Maxwell's 
equations is one of an incompressible elastic fluid. The 
obvious way to extend the validity of the model is to 
assume that the fluid is compressible. This means that 
it has to be a medium with vanishing compressibility, 
which happens when the dilational elasticity or viscos- 
ity coefficient is extremely large in comparison with the 
shear viscosity coefficient. As already mentioned, at 
this stage it is not clear whether the medium behaves 
as a liquid or solid when the compressional/dilational 
motions are considered. This question cannot be an- 
swered without staging an experiment in which the fab- 
ric of the metacontinuum could eventually be 'ripped' 
in order to settle this question. Actually, if the compres- 
sional/dilational motion is oscillatory, then it does not 
really make much difference if the metacontinuum is a 
solid or a fluid. As shown in the precedence, it has to 
be a fluid under shearing, but this does not impose any 
restriction on its behavior in the oscillatory compres- 
sion/dilation motions. Now, combining the two parts of 
the constitutive relation as given by Eqs. ( fTSI l the mo- 
mentum equation can be written as 

,dv „ ^ „ [(v H- 20divv Hquid, ^,„^ 
u(—+v-Wv) = -E + V ^' ^ (37) 

^dt \(A + 2T])divu solid, 

which is to replace Eq. ([TtI i. Eq. (ISTT l remains 
unchanged, while the incompressibility condition, 
Eq. (l22b . has to be replaced by the following: 



Tt 



-H V ■ Vyu = -//div V. 



(38) 



Alternatively, one can use the continuity equation in 
the fro m - J^^/j,- (see Chadwick (1999); Truesdell 
( Il965h ). where J is the determinant of the gradient of 
deformation tensor and fir is the constant density in 
the referential description. Note that, the form chosen 



above, is consistent with the Eulerian description. Note 
also that if the solid rheology is assumed in Eq. (iJTl l. 
then one has to add the defining equation for the veloc- 
ity components: v = du/dt, 

Eqs. ( ISTT i. ( l37b . and ( 1381 ) form a coupled system for 
the compressible metacontinuum, which is highly non- 
linear (even for the linear rheology) because of the de- 
pendence of the density ju on the motion. There are no 
conceptual difficulties to limit the model to the case of 
slight compressibility of the metacontinuum, and look 
for the propagation of harmonic compression waves 
therein. This raises the question about speed of the com- 
pressional waves ('sound') of the metacontinuum. In 
order to avoid ambiguous terminology we will not use 
the term sound when referring to these waves. Rather 
we will borrow a co i nage from the ancient school of 
Stoa (see lSamburskv ( 1965 )) and will call the compres- 
sional/dilational motions the pneuma. 

One obvious implication of the existence of waves of 
a different kind than the electromagnetic (shear) waves 
is that, in fact, there is more energy in the physical 
vacuum than that detected from electromagnetic inter- 
actions. Mechanically speaking, the pneuma waves 
are 'orthogonal' to the shear (EM) waves and they 
may not be detectable by devices based on electro- 
magnetic interactions. Indeed, they perfectly fit the 
bill of what is currently ca Ued 'dark energy' (see, e.g., 
Huterer and Tumeii (119991) ). One can begin thinking of 
how to detect pneuma waves only after some solutions 
for coupled compressional and shear waves of the sys- 
tem Eqs. (I2TI 1. dJTl l. and dSST l become available. 

Let us mention in the conclusion of this section 
that the above derived system is suitable for any 
(visco)elastic medium, such as th e ager-gelatin base 
phantom (see lCatheUne et al. (2003)). 



9. Conclusion: The Rational Continuum Mechanics 
of Electromagnetic Field 

The approach proposed here presents a self- 
consistent point of view based on the continuum me- 
chanics of the electromagnetic field whose shear de- 
formations are perceived as the phenomena of electro- 
magnetism. It is shown that the linearized governing 
equations of any incompressible elastic medium admit 
a 'Maxwell form'. Conversely, the Maxwell equations 
of electrodynamics are a strict mathematical corollary 
of the linearized governing equations of the incompress- 
ible Maxwell (visco)elastic fluid. 

This idea is further elaborated upon by deriving the 
frame indifferent formulation of the model of elastic flu- 
ids. It is argued that in some cases (as in the formu- 
lation of electromagnetism considered here and many 
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other technologically significant applications) the con- 
stitutive relation can be written for the stress vector 
rather than for the stress tensor, because what actually 
enters the momentum equations is the point-wise stress 
vector, which is the divergence of the stress tensor. 

From the frame indifferent governing system of elas- 
tic fluids, a 'Maxwell form' is derived that includes 
the terms of the original Maxwell equations along with 
terms stemming from the convective and convected in- 
variant time rates. These inertial terms are the progen- 
itors of the so-called laws of 'motional electrodynam- 
ics' : Biot-Savart, Oersted-Ampere, and Lorentz-force 
laws. The latter are usually assumed as additional em- 
pirical hypotheses to Maxwell's equations. This makes 
for a unified model of electromagnetism, which is truly 
covariant, by virtue of the fact that they are frame in- 
different. In other words, the new model is invariant to 
changes to frames that accelerate and deform. 

The continuum-mechanics formulation of the electro- 
magnetic field proposed here opens a new avenue of re- 
search cormected with the possible compressibiUty of 
elastic fluids. Consequently two kinds of Unear propa- 
gating waves can co-exist: shear waves (light, when in 
the visible spectrum) and compressional waves (called 
pneuma in this paper). It must be stressed that the speed 
of compressional waves is necessarily much larger than 
the speed of light. As a consequence, the metacontin- 
uum appears virtually incompressible to observes using 
tools designed to detect its shear motions (i.e, electro- 
magnetic phenomena). 
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